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ABSTRACT 

We present the spectrum of eigenfrequencies of axisymmetric acoustic- 
inertial oscillations of thin accretion disks for a Schwarzschild black hole 
modeled with a pseudo-potential. There are nine discrete frequencies, cor¬ 
responding to trapped modes. Eigenmodes with nine or more radial nodes 
in the inner disk belong to the continuum, whose frequency range starts 
somewhat below the maximum value of the radial epicyclic frequency. The 
results are derived under the assumption that the oscillatory motion is par¬ 
allel to the midplane of the disk. 

Keywords: Relativistic stars: black holes, structure stability and oscil¬ 
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1 ACOUSTIC-INERTIAL MODES 

We consider acoustic-inertial modes of oscillation in the inner part of an accretion 
disk, closely following the formalism of Nowak and Wagoner (1991, 1992). Trapping 
of the fundamental axisymmetric mode with no nodes in the vertical (z) direction 
was first demonstrated by Kato and Fukue (1980) in the Schwarzschild geometry. 
Nowak and Wagoner (1991) derive the equations of motion in a Lagrangian pseudo- 
Newtonian formalism and specialize to purely horizontal perturbed motions of the 
disk deriving eigenmodes and eigenfrequencies for the m = 0 (axisymmetric) and 
m = 2 (quadrupole) modes. Khanna et al. (2014) computed in an improved pseudo¬ 
potential the lowest radial modes (with up to three radial nodes) for azimuthal 
numbers m = 0 through to = 4. Here, we present the complete spectrum of 
horizontal axisymmetric acoustic-inertial disk modes in a pseudopotential which 
reproduces the properties of the Schwarzschild-metric epicyclic frequency (Kluzniak 
and Lee, 2002; Khanna et ah, 2014). The eigenfrequencies could be related to 
the quasi-coherent frequencies (QPOs) observed in the X-ray flux from black hole 
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and neutron star systems (for a review see van der Klis M., 2000), as well as in 
cataclysmic variables (Woudt and Warner, 2002, and references therein). 


2 EQUATION OF MOTION AND THE BOUNDARY CONDITION 

We model the Schwarzschild metric with a Newtonian pseudo-potential that repro¬ 
duces the Schwarzschild ratio of K^(r)/n^(r) = 1 — 6GM/{rc^): 

$kl(7’) = -(cV6) exp ■ (1) 

We have dropped an additive constant and renormalized the original Kluzniak and 
Lee (2002) potential by a factor of 1/e to guarantee the correct Schwarzschild value 
of n(r„is)- The orbital frequency can be obtained from n^(r) = r~^d^KL/dr, the 
radial epicyclic frequency from = (2n/r)d(r^Sl)/dr and the marginally stable 
orbit is at the zero of k, at r^s = QGM/{rc^). Fig. 1 compares our K^{r) with the 
Schwarzschild form and two other well-known pseudo-Newtonian models (Paczyhski 
and Wiita, 1980; Nowak and Wagoner, 1991). 

In this contribution we assume axisymmetric (m = 0) horizontal modes, with 
the perturbation vector in cylindrical coordinates (C*j?*,C*) = 0) exp(icrt). 



Figure 1. The Schwarzschild epicyclic frequency (squared) and its Newtonian models, 
from top to bottom: Paczyhski and Wiita (1980); Nowak and Wagoner (1991), our Eq. (1). 
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We use the equation of motion for 'l'(r) = \/^fPrff{r) derived in the Lagrangian 
formalism of Friedman and Schutz (1978) by Nowak and Wagoner (1991), 

Cgd^5'/dr^ + (cr^ - k^)^' = 0, 


who also show that in the WKB approximation the azimuthal component of the 
equation of perturbed motion for thin disks reduces to = 2i{p./a)ff. 

Following Khanna et al. (2014) we rewrite the equation of motion, and the bound¬ 
ary condition that the Lagrangian perturbation of pressure vanishes at the unper¬ 
turbed boundary, in dimensionless form as 



with the boundary condition at a; = 0 


d^* _ 4' 

dx 2 ’ 


( 3 ) 


where H is the half-thickness of the disk, a = rms, and the dimensionless variables 
are given by r = a(l + x), a = cr/n(a), k(x) = K(r)/Vl{a). The speed of sound 
Cs = VtF/p was eliminated with the condition of vertical hydrostatic equilibrium. 
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Figure 2. Two radial overtones for m = 0 horizontal oscillations of a thin {H/a = 10”®) 
accretion disk. Left Panel: A trapped oscillation with p = 8 radial nodes. Right Panel: 
An oscillation penetrating the epicyclic barrier (with fi = ^ radial nodes in the inner 
accretion disk and an unlimited number of radial nodes in the outer disk). Plotted are the 
radial wavefunction oc : solid blue line (left scale); eigenfrequency (squared) and 

the epicyclic frequency (squared), both normalized to orbital frequency at the inner edge 
of the disk, (j^/fl^(rms): dashed-dotted green line, K^(r)/fl^(rms): dashed red line (right 
scale). 
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3 THE EIGENFREQUENCY SPECTRUM 

We have numerically solved the eigenvalue problem given by Eqs. (2), (3), for a 
thin disk of H/a = 10“^, and present in Table 1 the eigenfrequencies for modes 
with /i = 0,1,9 radial nodes in the inner disk. The lowest nine eigenfrequencies 
(/X = 0 through 8), exhausting the discrete spectrum, correspond to oscillations 
which are trapped in the inner disk. As already noted by Kato and Fukue (1980), 
for eigenfrequencies exceeding the maximum of the epicyclic frequency, > K^ax) 
the acoustic wave ranges throughout the disk (see also Kato et ah, 1998), these 
frequencies belong to the continuum spectrum. The tenth entry in Table 1, with 
/X = 9 radial nodes in the inner disk, also belongs to the continuum, although it has 
a frequency below the maximum of the epicyclic frequency (Figs. 2, 3). 

In Fig. 2 we present the two eigenmodes corresponding to the last two entries in 
Table 1. The equations being linear in dx, we normalize the wavefunction to unity 
at the inner edge of the disk: 4'(rms) = 1 for illustration purposes. The left panel 
shows the highest-frequency eigenmode in the discrete portion of the spectrum of 
axisymmetric (m = 0) horizontal disk oscillations, the wavefunction of this mode 
has /X = 8 radial nodes. Note that the wave becomes evanescent for thus 

trapping the /x = 8 mode to the left of Kmax- Some of the lower overtones have been 
illustrated in Khanna et al. (2014). 

The right panel of Fig. 2 illustrates one of the lowest frequency modes in the 
continuum. Here, a < Kmax and is so close in value to Kmax that the wave is 
transmitted through the epicyclic barrier to the outer disk, where it has an unlimited 
number of radial nodes in addition to the /x = 9 radial nodes in the inner disk. As 
far as we are aware, this is a new finding, which has never been reported before. It 
may have an interesting astrophysical consequence. If the oscillations arise close to 
the marginally stable orbit, as suggested by Paczyhski (1987), the ones transmitted 
to the outer disk are likely to be more easily observable, in that they may modulate 
the emission from large parts of the disk. 

Table 1: Spectrum of eigenfrequencies for horizontal modes 


mode 

m 

radial nodes 

eigenfrequency 

a X \/(a^/GM) 

mode 

status 

ref. 

0 

0 

0.098829oooo 

trapped 

Khanna et al. (2014) 

0 

1 

0.172137256O 

trapped 

Khanna et al. (2014) 

0 

2 

0.2049091375 

trapped 

Khanna et al. (2014) 

0 

3 

0.22691222OO 

trapped 

Khanna et al. (2014) 

0 

4 

0.2433069ooo 

trapped 

this work 

0 

5 

0.2561130000 

trapped 

this work 

0 

6 

0.2663460000 

trapped 

this work 

0 

7 

0.2745830000 

trapped 

this work 

0 

8 

0.2811589OOO 

trapped 

this work 

0 

9 

0.2862229OOO 

not trapped 

this work 
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4 DISCUSSION 

We consider accretion disk oscillations in a Newtonian model of the Schwarzschild 
metric, Eq. (1), which accurately models the radial epicyclic frequency, at least 
close to the marginally stable orbit, see Fig. 1. No model is perfect, so although 
we correctly reproduce the ratio of epicyclic to orbital frequency K(r)/n(r) = 
y/T^^'6GMjJr^, and the correct value of orbital frequency at the marginally 
stable orbit, r2(rms) = c^/(-\/216 GM), the maximum of k occurs at r = (3 + 
■\/2 T)(GM/c^) ~ 7.58(GM/c^) instead of the Schwarzschild value r = 8{GM/c^). 
Further, the equations of motion for the oscillation of the disk fluid were derived in 
a Newtonian formalism, not in full GR. These departures from GR may limit the 
quantitative accuracy of the presented results when applied to real black hole (or 
neutron star) accretion disks. An additional assumption which may not be quite 
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Figure 3. The spectrum of m = 0 (axially symmetric) horizontal oscillations of a thin 
{H/a = 10“®) accretion disk for the potential of Eq. (1). Plotted are the epicyclic fre¬ 
quency (dashed red curve), and the eigenfrequencies a in the discrete set (dashed-dotted 
green lines) and in the continuum (shaded blue region). All frequencies were scaled with 
[GM. Here, and throughout the paper, a = rms. 
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accurate is that the oscillations of the disk are strictly parallel to the midplane of 
the disk, i.e., that the perturbation vector has a zero vertical component, = 0. 

We find that the spectrum of horizontal oscillations is composed of nine discrete 
frequencies and a continuum (Fig. 3). For the discrete spectrum the wave propa¬ 
gation region corresponds to those regions where > K^{r) and is separated into 
the inner region of trapped oscillations, from r = rms to r « (7/6) rms, and an 
outer region extending to r >> rms (Kato et ah, 1998). However, the lowest fre¬ 
quency modes in the continuum, which satisfy cr < Kmax, are transmitted through 
the epicyclic barrier, and thus fluctuations in the inner disk may be transmitted to 
the outer disk for frequencies close to the maximum of the epicyclic one, « «^max- 
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